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Abstract 

The  behavior  of  multi-dimensional  discrete  Boltzmann  systems  with  highly  oscillatory  data  is 
studied.  Homogenized  equations  for  the  mean  solutions  are  obtained.  Uniform  convergence  of  the 
oscillatory  solutions  of  the  discrete  Boltzmann  equations  to  the  solutions  of  the  corresponding 
homogenized  equations  is  established.  Moreover,  we  find  that  the  weak  limits  of  the  oscillatory 
solutions  for  a  model  of  Broadwell  type  are  not  continuous  functions  of  the  discrete  velocities. 
Generalization  of  the  above  results  to  problems  with  multiple-scale  initial  data  is  also  established. 

X 

•A, 

1.  Introduction 


One  of  the  essential  problems  in  nonlinear  partial  differential  equations  is  to 
understand  how  the  behavior  in  the  microscopic  level  affects  the  behavior  in  the 
macroscopic  level.  One  approach  is  to  find  the  corresponding  effective  or  ho¬ 
mogenized  equations  which  determine  the  weak  limits  of  the  oscillatory  solutions 
(see  [3]),  Such  homogenization  results  are  important  both  for  the  understanding 
of  the  nonlinear  interaction  between  the  high  frequencies  and  for  the  study  of  the 
numerical  approximation  for  problems  with  oscillatory  solutions  (see  [7],  [8],  [9]). 

We  choose  to  study  the  nonlinear  discrete  Boltzmann  equations  in  kinetic 
theory  of  discrete  velocity.  In  such  equations,  high  frequency  components  can  be 
transformed  into  lower  frequencies  through  nonlinear  interaction,  thus  affecting 
the  average  of  solutions.  In  this  paper,  we  study  the  homogenization  theory  of  the 
discrete  Boltzmann  equations  in  multi-dimensional  space  and  with  finitely  many 
real-valued  velocities  (see  equations  (4.1)).  We  assume  that  the  initial  values  are 
of  the  form  a(x ,  x/e)  with  a(x,  y)  l-periodic  in  each  component  of  y.  Our 
results  show  that  the  behavior  of  oscillatory  solutions  for  a  model  of  Broadwell 
type  (see  equations  (2.1))  is  very  sensitive  to  the  velocity  coefficients.  It  depends 
on  whether  a  certain  ratio  among  the  velocity  components  is  a  rational  number  or 
an  irrational  number.  Furthermore  we  find  that  the  form  of  homogenized 
equations  depends  on  the  velocity  coefficients,  and  the  weak  limits  of  the 
oscillatory  solutions  are  not  continuous  in  the  velocity  components.  This  singular 
behavior  for  a  model  of  Broadwell  type  is  not  shared  by  the  simple  Carleman 
model  (see  equations  (3.1)). 

Our  study  also  shows  that  the  structure  of  oscillatory  solutions  for  a  model  of 
Broadwell  type  is  quite  stable  (in  the  sense  of  Theorem  2.3)  when  we  perturb 
velocity  coefficients  around  irrational  numbers.  In  this  case,  the  resonance  effect 
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of  u  and  v  on  w  vanishes  in  the  limit  of  e  -*  0.  However,  the  behavior  of 
oscillatory  solutions  for  a  model  of  Broadwell  type  becomes  singular  when 
perturbing  around  integer  velocity  coefficients.  There  is  a  strong  interaction 
between  the  high  frequency  components  of  u  and  v,  and  the  interaction  in  the  u, 
v  term  would  create  oscillation  of  order  0(1)  on  the  w  component  even  in  the 
limit  as  £  -*  0.  In  [14],  Tartar  showed  that  for  the  Carleman  model  the  weak 
limits  of  all  powers  of  the  initial  oscillatory  data  will  uniquely  determine  the 
weak  limits  of  the  oscillatory  solutions  at  later  time.  We  found  that  this  is  no 
longer  true  for  a  model  of  Broadwell  type  with  integer-valued  velocity  coeffi¬ 
cients. 

The  homogenization  theory  of  the  Carleman  and  Broadwell  models  with 
oscillatory  initial  data  has  been  studied  by  McLaughlin,  Papanicolaou  and  Tartar 
[11].  They  proved  that  the  oscillatory  solutions  of  the  Carleman  and  Broadwell 
models  converge  strongly  in  Lp-norm,  p  ^  oc,  to  the  solutions  of  the  corre¬ 
sponding  homogenized  equations.  By  using  certain  ergodicity  property  of  the 
oscillatory  solutions  and  taking  into  account  cancellations  among  high  frequency 
components,  we  are  able  to  obtain  homogenization  results  for  more  general 
discrete  Boltzmann  equations.  Moreover  we  establish  uniform  convergence  of  the 
oscillatory  solutions  to  the  solutions  of  the  corresponding  homogenized  equa¬ 
tions.  This  uniform  convergence  result  is  essential  in  the  convergence  analysis  of 
particle  methods  for  the  discrete  Boltzmann  equations  (see  [7],  [8],  [9]). 

The  paper  is  organized  as  follows.  In  Section  2,  we  study  a  model  of 
Broadwell  type  in  detail  and  compare  the  homogenization  results  with  those  for 
the  Carleman  model.  Section  3  contains  homogenization  results  for  the  problems 
in  which  the  initial  data  are  of  more  than  two  scales.  In  Section  4,  we  extend  the 
results  of  Section  2  to  the  discrete  Boltzmann  equations  in  multi-dimensional 
space  with  finitely  many  velocities. 

2.  Behavior  of  Oscillatory  Solutions  in  the  Model  of  Broadwell  Type 

The  Broadwell  model  describes  a  three-dimensional  model  of  rarefied  gas  in 
which  particles  travel  with  speed  c  in  either  direction  along  a  coordinate  axis  (see 
[4]).  If  particles  traveling  in  opposite  directions  collide,  they  are  equally  likely  to 
move  in  each  of  the  three  coordinate  directions  after  collision,  with  velocities  of 
opposite  sign.  Other  collisions  can  lead  to  an  exchange  of  velocities.  We  denote 
by  N+(x,  y,  z,  t)  the  number  density  of  particles  with  velocity  (c,0,0);  a  similar 
notation  is  used  for  N2±  and  JVj*.  Then  the  resulting  equations  are 

~  c~Jx~  =  5*(AW+  2AW), 

etc.,  where  o  is  the  frequency  of  collision. 
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Here  we  consider  the  special  case  of  one-dimensional  motions  in  which  the 
TV’s  are  independent  of  y,  z,  and  furthermore  TV2+  =  N{  =  7V3+  =  N{.  Setting 
Ar1+=  u(x,  l),  Ni=v(xyt),  TV2+  =  w(x,  t)  and  rescaling  the  variables  so  that 
c  —  1,  a  =  §,  we  then  obtain  the  1-D  Broadwell  equations 


du 

du 

9  r\ 

dt 

+  J^  +  UV- 

w2  =  0, 

dv 

dv 

Tt 

-Ji+uv~ 

w2  =  0, 

dw  ,  ■>  a 

-jp —  uv  +  w  =  0. 


A  lot  of  effort  has  been  made  in  obtaining  global  solutions  for  the  discrete 
Boltzmann  equations.  For  the  1-D  Broadwell  model,  Nishida  and  Mimura  [12] 
first  showed  that  a  global  solution  exists  when  the  initial  values  are  small  in  some 
sense.  Their  result  has  been  generalized  to  more  general  1-D  discrete  Boltzmann 
equations  by  Crandali  and  Tartar  [13],  Cabannes  [5],  lllner  [10]  and  Beale  [2], 
among  others. 

To  study  how  oscillatory  solutions  depend  on  velocity  coefficients,  we  intro- 


duce  an  additional  term 
model.  We  get 

a  dw/dx,  jaj  <  1,  in  the  last  equation  of  the  Broadwell 

(2.1a) 

du  du  2 

17 +  Tx+uv~w  =  °’ 

(2.1b) 

dvdv  2 

=  o. 

(2.1c) 

dw  dw  2 

17  +  a17~uv  +  w  =0- 

We  refer  to  equations  (2.1a-c)  as  a  model  of  Broadwell  type.  Suppose  in 
addition  that  initial  values  are  given  by 

(2.2)  u(x,0)  =  u0(*»  j),  u(x,0)  =  «0(x,  y),  w(x,0)  =  j). 


where  we  assume  that  u0(x,  y),  v0(x,  >>)  and  y)  are  1-periodic  in  y. 

We  denote  by  ue,  v(  and  wt  the  solutions  of  equations  (2.1)  and  (2.2). 

The  behavior  of  solutions  ue,  vc  and  *vt  as  e  -»  0  is  very  sensitive  to  the 
coefficient  a.  This  is  described  by  the  following  homogenization  result. 
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Case  I:  a 

=  m/n  (m  and  n 

(2.1a-c)  are 

(2.3a) 

du  dfJ 

dt  +  dx 

(2.3b) 

dV  dV 

1  x 

1  |<^> 

1 

-/>>=<*• 


(2.3c) 


SW  3W 

sr*^*w 


-\£U[X'  y  +  (?  "  T)Z’  'M*’  '  +  (7  +  ^  ')  =  °‘ 

If  m  =  0,  then  n  =  1  in  (2.3c). 

Case  II:  a  is  an  irrational  number;  the  homogenized  equations  become 


(2.4a) 

(2.4b) 


dU 

dt  ^ 

dx 

dv 

dv 

dt  " 

dx 

(2.4c)  ^  +  W2  ~  (jjV*.  J'.  0  *  0  =  °- 

The  initial  data  in  both  cases  are  given  by 

U(x,  y, 0)  =  u0(x,  y),  V(x,  y,0)  =  v0(x ,  y), 

IT(x,  y,0)  =  w0(x,  y). 

Here  we  have  assumed  that  smooth  and  bounded  global  solutions  of  (2.3a-c) 
and  (2.4a-c)  exist  up  to  time  T. 


(2.5) 


Theorem  2.1.  For  smooth  and  bounded  non-negative  initial  data .  the  solutions 
of  (2.1a-c)  and  (2.2)  converge  to  those  of  the  corresponding  homogenized  equations 


strongly  in  the  L^-norm, 

ue(x,t)  l/(x,  e  ,tj 

-  0 

as 

e  -*  0, 

ve(x,t)  V^x,  £  ,tj 

-  0 

as 

£  -*  0, 

*e(x’')  W(X’  e  ’*) 

-*  0 

as 

e  -»  0  for  0  S  (  g  r. 
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Remark  2.1.  The  homogenization  result  above  can  be  generalized  to  the  case 
when  uQ(x,  y),  v0(x,  y)  and  w0(x,  y)  are  periodic  functions  in  y  with  arbitrary 
periods. 

Remark  2.2.  In  the  case  when  a  =  0,  our  homogenization  result  is  the  same 
as  that  obtained  in  [11].  The  only  difference  is  that  we  obtain  uniform  conver¬ 
gence  instead  of  convergence  in  L^-norm,  p  <  oo,  of  [11]. 


Remark  2.3.  The  local  existence  result  of  the  homogenized  equations  (2.3)  or 
(2.4)  can  be  obtained  by  classical  analysis  for  smooth  and  bounded  non-negative 
initial  data.  The  global  existence  result  for  the  homogenized  equations  then 
follow  by  combining  the  known  global  existence  results  for  (2.1)— (2.2)  (e.g.  [10]) 
with  Theorem  2.1.  Therefore  the  value  of  T  in  Theorem  2.1  is  arbitrarily  large. 


Lemma  2.1.  Let  /(x),  g(x,  y)  e  C1.  Assume  further  that  g(x,  y)  is  n-peri- 
odic  in  y  and  satisfies  the  relation  fog(x,  y )  dy  =  0.  Then  for  any  constants  a  and  b, 
we  have 


1 

b  -  a 


g  Ce. 


Proof:  Express  g(x,  x/e)  as 

(2.6)  g(x,  j)  =  7)  ds  -  f )  ds. 

Since,  for  any  real  number  d, 

C"j f (*’  y)  y)dys0> 

we  conclude  that 

\fj(X'  l)  *  Q£’  e)  *|  =  C2£- 

From  this,  we  deduce  that 

■=7  £/<*)»(*.  t)*| 

-  — 1 — |/(6)^k(*.  ;)*  -  jf*(jf  *(*. — 


I  b  ~  a  | 


dx 


dx 


^  Ce. 


This  completes  the  proof  of  the  lemma. 
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Lemma  2.2.  Suppose  that  f(x,  y,  z)  is  continuous.  Moreover ,  we  assume  that / 
is  \-periodic  in  y  and  z.  If  is  an  irrational  number ,  then  we  have 


vh £('(*■ ^)- -/;/>■  *'>♦*)  - 

as  e  ->  0  for  any  constants  a,  b,  a,  and  x2. 


^  C(e)  -  0. 


Proof:  We  assume  for  simplicity  that  /  has  been  normalized  so  that 


We  first  prove  the  lemma  in  the  case  when  f(x,  y,  z)  is  independent  of  the  first 
variable  x.  By  change  of  variables,  we  can  further  reduce  the  problem  to  showing 
that 


1  fT 

Y  J  f(x,xl2  +  Xx)dx 


0  as  T  -*  oo. 


where  A  =  y2/y,  is  an  irrational  number,  T  —  yfb  -  a)/e,  xh2  =  (x2  -  Xx^/e. 
We  need  only  to  show  this  for  integer-valued  T.  Note  that 


fN/(x,xu2  +  \x)  dx  =  jf  £  (n  +  lf{x.x12  +  Xx)dx 
0  0  n 

(2.7) 

1  /v_1  l 

=  -77  Y,  l  f(x ’  a,  2  +  «A  +  Ax)  dx. 

^  n- 0J0 

Define 

Hy)  =  jfVu*  y  +  dx- 


Then  F(y)  is  a  1-periodic  continuous  function.  Applying  the  well-known  theo¬ 
rem  of  equipartition  modulo  1  (Bohl-Serpinskii-Weyl)  of  ergodic  theory  (see  [1]), 
we  obtain 


(2.8) 


N-  1 


Tv  £  ^(*1,2  +  «A)  -  pF(y)  dy 


n  —  0 


0  as  N  oo. 


independently  of  the  value  of  xl  2.  By  interchanging  the  order  of  integration,  we 
get 


PF(y)  dy  =  ppfiy* z)  dydz  =  °- 
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Thus  (2.7)  and  (2.8)  yield 


1  M  A  ■*!  +  Yi*  *2  +  ' 


hiu 


dx  -»  0  as  e  -»  0. 


For  the  case  where  /  depends  on  its  first  variable,  we  have,  in  light  of  (2.6), 


(2.10) 


Note  that  <9/(*,  >>,  z)/dx  is  1-periodic  in  the  j  and  2  variables,  and  satisfies 


y> z)  dydz  =  ^/070l/(x’  >'■ 2)  dydz  =  0> 

by  the  assumption  on  f(x,  y,  2).  Applying  (2.9)  to  the  integrals  with  respect  to  ds 
on  the  right-hand  side  of  (2.10),  we  obtain 

This  completes  the  proof  of  Lemma  2.2. 

Proof  of  Theorem  2.1:  We  only  give  the  proof  for  the  case  when  a  is  an 
irrational  number.  The  case  when  a  =  m/n  follows  similarly  by  using  Lemma 
2.1.  Subtracting  equation  (2.4c)  from  equation  (2.1c)  and  integrating  the  resulting 
equation  along  their  characteristics  from  0  to  t ,  we  obtain 


wt(x,t)  -  H'j*,  * 

=  —  J  ^wt(x  —  a(t  —  s),  s)2  —  W^x  -  a(t  —  s).  ~~~~~ 
(2.11)  +  T(x .  0  +  -  a(t  -  s ),  -  «(/  -  j),  r) 


r  dr 


-  J'i/lx  -  «(/  -  s). 


(a  -  at)  +  (a  —  l)s 


-  a(t  “  *), 


(x  ~  at)  +  (a  y  l)s 


,  s|  ds. 
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where  T(x,  t)  is  defined  by 


T(x,  t )  =  jff/|x  -  a(t  -  s), 


x  -  at  (a  -  l)s 

-  +  - - —  .  5 

e  e 


,  x  -  at  (a+  l)j  \ 

•F^x-«(l-  s),— —  +  i s  j  ds 

-  JoIJqu(x  ~  a(i  - s )’ z’ s )  -  •y)>  ■ 

For  fixed  x  and  t,  define  A  *  (a  +  l)/(a  -  1)  and 

f(s,  y,  z)  =  f/(x  -  a(/  -  s),  y,  s)V(x  -  a(t  -  s),  z,  s), 
*1  =  x2  =  x  -  at,  Yi  =  a  -  1,  y2  =  a  +  1. 

Then  the  T(x,  t )  term  in  (2.11)  becomes 

/;(/(»■ 

which  is  bounded  uniformly  by  C(e)  by  Lemma  2.2. 

Define 

A/=  sup  f  sup  {|«e(x,  01,  |t/(x,  y,  01,-  •  |»Xx,  y,  t)|}\ 

t  v  *.  y  ' 

We  deduce  from  (2.11)  that 


s )  dy  ds 


(2.12) 


ve(x,  /  )  -  W^x,  —  -  at ,  t  j  S  2 M  J  g{x  -  a(t  -  s),  s)  ds  +  C(e), 


where  G(x,t)  is  defined  by 


G(x,  t)  =  ut(x,  t )  -  f/|x,  +  °^x' r)  _  V{x'  ^7“^’ f) 

+  0  -  w[x,  X  -  f  j  . 


(2.13) 
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Similarly,  we  obtain 


(2.14) 

(2.15) 


ue(x,  t)  -  (/|jc,  —  j —  ,/j  ^ 2M  j  G{x  -  t  +  s,  s)  ds  +  C{e ). 


vt(x,  t)  -  V^x,  —  -  - ,  r  j  ^2M  j  G{x  +  t  -  s,  s)  ds  +  C{e). 

Adding  (2.12),  (2.14)  and  (2.15)  yields 

G{x,  t)  ^  2M J  (G(x  -  t  +  s,  s) 

(216)  +  G(x  +  t~s,s) 

+  G(x  -  a(t  -  s),  s))  ds  +  3C(r). 


Define 


£(/)  =  sup  { G (  a ,  / ) ) . 


It  follows  immediately  from  (2.16)  that 

(2.17)  £(/)  g  6M^E(s)ds  +  3C(e). 


Application  of  the  Gronwall  inequality  to  (2.17)  then  proves  the  theorem. 

Suppose  that  w0  is  y  independent.  The  oscillation  of  u  and  i-  will  create 
oscillations  on  w  at  later  time.  The  homogenized  equation  (2.3c)  indicates  that  w 
remains  oscillatory  as  e  -»  0  if  a  is  rational.  However,  if  a  is  an  irrational 
number,  (2.4c)  implies  that  W(x,y,t)  is  y  independent.  Thus  we  expect  that 
there  is  some  kind  of  singularity  in  the  high-order  powers  of  solutions.  Since  the 
equations  are  nonlinear,  such  a  singularity  would  affect  the  local  average  of 
solutions.  This  is  described  by  the  following  theorem. 


Theorem  2.2.  Let  a0  =  m/n.  Assume  w(x,0)  =  w0(x)  and 

y  +  {j~  y +  (?  +  dz  *  °- 

Then  at  least  one  of  the  following  limits  does  not  hold  as  a  -»  a0: 

fva  dy  —  PG  dy,  f  K  dy  -  f  V  dy .  f  Wa  dy  -  f  W  dv . 

J  0  J0  J0  J0  J0 


where  Ua ,  Va  and  Wa  are  solutions  of  equations  (2.3)  or  (2.4). 
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Proof:  It  can  be  shown  that 


£»l"u*(x’  ,*'u~  ')*•  '  K(*-  y + (7 +  ')2- ') <feA 


-  uv)uv)- 

Let  U{m)  =  jQUmdy.  For  irrational  a,  we  first  integrate  equation  (2.3c)  corre¬ 
sponding  to  a0  and  equation  (2.4c)  corresponding  to  a  in  y  from  0  to  1  and  then 
integrate  from  0  to  t  along  their  characteristics,  respectively.  The  difference  of  the 
resulting  equations  gives 

W<l)(x,t)  -  tVa(I>(x,t)  +  f\w0(x  -  a0t ,  y)  -  w0(x  -  at,  y))  dy 

Jo 

(2.18)  + 

=  -  ~  a(>  ~  s).s)  -  H^>(x  -  a0(t  -  s),s))  ds. 


Since  y,  0)  is  ^-independent,  ITa(x,  t )  is  independent  of  y  for  irrational 

a  according  to  equation  (2.4c).  Thus  WJ2)  =  (lTaa>)2.  On  the  other  hand,  the 
assumption  on  u0,  v0  and  equation  (2.3c)  imply  that  IF  ( x,  >\  O  is  a  nontrivial 
function  of  y  at  least  in  some  interval  (0,  fj  with  f,  >  0.  Thus  we  have 


(2.19) 


j[V.“  -  <')  <*  -  f0"((^"f  -  K")') * 

'K"-Kn)!K 


Since  Wag( x ,  y,  t)  is  not  independent  of  y,  Schwartz’s  inequality  yields 

K'  ~  «'  f  >  «• 

Therefore  the  last  term  on  the  right  of  (2.19)  is  a  non-zero  function  independent 
of  a.  Consequently, 

f4<l)  -  (C  v±Vl  -  C’  w«l)  -  K]  as  «  -  «0’ 

would  contradict  (2.18)  and  (2.19).  This  completes  the  proof  of  Theorem  2.2. 


The  situation  for  a0  irrational  is  quite  different.  Weinan  E  and  I  can  show 
that  the  solution  is  ‘structurally  stable’  near  a0  in  the  following  sense. 
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Theorem  2.3.  Suppose  that  solutions  Lla(x,  \  t)  and  Va(x,  *,  /)  of  (2.3)  or 
(2.4)  belong  to  C 2.  If  a0  is  an  irrational  number ,  then  solutions  Ua.  Va  and  Wa  of 
(2.3)  or  (2.4)  are  continuous  functions  of  a  at  a0.  Moreover,  for  any  positive 
integer  m,  we  have 

lim  flUa(x.  y,t)m  dy  =  C U  (x,  y,t)m  dy, 

lim  f'vjx,  y,  t)m  dy  =  fv '  (x.  y,  t)m  dy, 

a-a0J  0  J0 

lim  fl  IVJx,  y,t)m  dy  =  f  y,t)m  dy. 

a  —  a0J  0  •'() 

Proof:  Case  1:  a  a0  and  a  =  m/n.  We  assume  that  0  <  |o0|  <  1.  and 
that  m  and  n  are  mutually  •’rime  with  n  >  0.  Since  a0  is  an  irrational  number, 
clearly  we  have  lima_ao«(a)  =  oo. 

Define  the  function  H(x,  y ,  t)  by 

s  \[u\x'y  +  (?  -  >■  +  (7  +  \)z,t)dz 

~(/0>u°(x'  *0*)(jjV.<*.  x>()dy )■ 

We  first  show  that  limn  JXf/(x,  y,t)  =  0.  Since  U„(x.  y.  t)  and  Kn(x.  y,  r  1  are 
1 -periodic  functions  in  y,  we  can  expand  Va  and  Va  by  their  Fourier  series.  We 
get 


(2.20) 


H  =  \  f  (  L  a*exp{2irifc(  y  +  (m/n  -  l)r )}) 
HJo'k+o  ' 


fi/exp{27r//(  y  +  ( m/n  +  l)r)}  dz 


It  can  be  shown  that  dfUa  is  bounded  independently  of  a.  Therefore,  we  have 
\ak\  ^  c(\/k2)  and 

53  a kexp{2nik(  y  +  ( m/n  -  l)c)}  <1  C’/n. 

|A|>n/3 

Similarly,  we  have 

5!  fi,exp{27r//(  v  +  (m/n  +  1)j)}  ^  C/n. 

I  h>n/i 

On  the  other  hand,  note  that,  for  0  <  |/cj  <,  \n,  0  <  |/|  ^  \n. 


(2.21) 


k(m/n  —  1)  +  l(m/n  +  l)  *  0. 
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Suppose  otherwise.  We  then  get  m/n  =  (k  -  l)/(k  +  /).  But  \k  +  /|  <;  \n  <  n , 
| m\  <  n,  which  contradicts  the  assumption  that  m  'nd  n  are  mutually  pnme. 

As  a  result  of  (2.21),  we  can  easily  show  by  interchanging  the  order  of 
integration  and  summation  that 


\i  {  E  a*exp{2fl7*(  y  +  (m/n  -  l)r)} 

0  ^0<|*|s«/3 


E  b,tx^{2vil(y  +  (m/n  +  l)z)} 
0<|/lS«/3 


=  0. 


Therefore  we  have  proven  that  Iim„ y,  t)  =  0  bv  showing  that 
(2.22)  \H(x.y,t)\gC/n, 


where  C  is  independent  of  a. 

Integrate  equation  (2.4c)  corresponding  to  a0  and  equation  (2.3c)  correspond¬ 
ing  to  a  from  0  to  t  along  their  characteristics,  respectively.  The  difference  of  the 
resulting  equations  gives 

Wa(x,  v.t)  -  Wjx,  y.f)  =  «0(x  -  at,  y)  -  wa(x  -  a()t.  y) 

+  j'va(x  -  a(t  -  s),  y.s)  dy  )|  jf’  Va( x  -  <i(  t  -  s ).  y.  5 )  dy  j  ds 

(2.23)  -  jT'(/oV*  -  a»('  ~  s )'  >'  ? )  dv  j ( £ !"„„( X  -  «„(/  -  s),  y,  s)  dy)  ds 

+  (  H(x  -  a{ t  -  s ),  v,  s )  ds 

M) 

+  ['{KJ.X  ~  ao(<  -  s).y.sY  ~  WJx  -  a(t  -  s).  y.  r)‘)  ds. 

Jo v 

Define 


G(x,t)  s  ||t/o(x,  \{)  -  U  (x,  *,  /)||t*  +  HK,(  x,  •.  t)  -  T  (x,  *,  t)|i/> 


+  ll^«(x, -,t)  -  ^(*.*,r)||#-. 
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Since  Ua,  Va  and  Wa  are  bounded  and  u0,  v0  and  w0  have  bounded  partial 
derivatives,  we  can  show  that 


(2.24) 


\Wa(x,  y,  t)  -  W  ( jc,  y,  f)|  ^  2 A/  (‘g(x  -  a(t  -  s),t)  ds 
>  •'o 


+  C/n  +  C|a  —  a0|. 


Similarly,  we  have 

| U„(x,y.t)  -  Uan(x,  y,  f)| 

^  2A { j  G(x  —  (t  —  5).  /)  ds  +  C/n  +  C|a  -  a0|, 

\Kix-  .v.O  -  KS*.  v.0| 

^  2M  j  G(x  +  (t  -  s ),  t)  ds  +  C/n  +  C|a  -  a0|. 


(2.25) 


(2.26) 


Define  E(t)  =  sup.GU,  1).  Adding  (2.24),  (2.25)  and  (2.26)  yields 

(2.27)  E(t)  s  6 A/  ('e(s)  ds  +  3C/n  +  3C]a  -  a0|. 

•'0 

The  Gronwall  inequality  then  implies  the  theorem  for  Case  1. 

Case  2:  a  is  an  irrational  number.  In  this  case  the  homogenized  equations 
corresponding  to  a  and  a0  are  of  the  same  form  (2.4).  The  proof  is  identical  to 
the  second  step  in  the  Case  1  beginning  from  (2.22).  Hence  the  proof  of  Theorem 
2.3  is  complete. 


The  proof  of  Theorem  2.3  contains  the  following  result. 


Theorem  2.4.  If  a„  =  m/n  and  n  »  1.  then  there  exists  5(n)  >  0  such  that, 
for  |a  -  a0|  <  S(n), 


f'ujx,  y,  t)k  dy  -  flU  (x,  y,t)k  dy 
•'ll  Jo 


+A(a 


Similar  expressions  hold  for  V  and  W.  Here  A(8)  -»  0  as  8  -*  0.  and  k  is  any 
positive  integer. 


The  following  theorem  tells  us  how  the  structure  of  oscillatory  initial  da  j 
affects  the  higher-order  powers  of  solutions  in  the  future.  For  the  Carleman 
model,  Tartar  [14]  has  shown  that  the  weak  limits  of  u/x.  t)m  are  uniquely 
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determined  by  the  weak  limits  ot  uf(x,0)k  for  all  k.  For  a  model  of  Broad  well 
type,  we  show  that  this  is  true  only  if  «  is  an  irrational  number. 


Theorem  2.5.  Case  I:  a  is  an  irrational  number.  Assume  {Uu,  Va,  Wu)  and 
{Uh,  Vh,  Wh }  are  two  different  solutions  of  system  (2.4).  Then  we  have 

(2.28a)  C  Uu(x,  y,  t)m  dv  =  C Uh(x.  y.  t)m  dy. 

■'ll  J(> 

(2.28b)  f  ‘  Va(x,  v,  t)m  dy  =  Vh(x.  v. /)"  dy. 

Jo  Jo 

(2.28c)  flfVa(x,  y,  t  )m  dy  =  C  Wh(x.  y.  t)m  dy. 

for  all  integer  m  ^  1  and  0  <  t  jg  T provided  that  (2.28a).  (2.28b)  and  (2.28c)  are 
valid  at  t  =  0  for  all  integer  m  ^  1 . 

Case  II:  a  is  a  rational  number.  Let  {Ua.  Va.  Wa }  and  { Uh.  Vh.  Wh }  be  two 
different  solutions  of  system  (2  3).  Then  (2.28a).  (2.28b)  and  (2.28c)  may  not  be 
true  in  general  for  all  integer  m  S  1  and  0  <  t  g  T  even  if  (2.28a).  (2.28b)  and 
(2.28c)  are  valid  at  t  =  0  for  all  integer  m  ^  1. 


Proof:  Case  I:  We  follow  closely  the  proof  given  by  Tartar  in  [14]. 
Throughout  the  proof  we  shall  use  the  notation  U(m)(x.  t )  =  f^U(x.  y.  t)m  dy. 
For  irrational  a.  one  could  easily  derive  equations  for  L'J m \  Vjm>  and  IVJ '  as 
follows: 


(2.29a) 

(2.29b) 

(2.29c) 


dW‘ 


dt 


dU'm) 

dU'm) 

mUf' 

dt 

+  dx 

+ 

mUJm)Vju  - 

"1TJ-' 

=  0. 

dV(m' 

dv{m) 

mVJm>UJl)  - 

a 

dt 

_ _ a 

dx 

+ 

mvr 

Du/(2) 

a 

=  0. 

m) 

dW'm) 

Wf"'u  -  ml[ 

1  > 
a 

—  +  a 

“  4- 

dx 

m 

lOl’ll) 

j  a 

=  0: 

U*m>,  Vlm)  and  Wlm)  satisfy  the  same  equations  (2.29a-c). 

From  the  global  existence  results  [12J.  [13]  and  [2],  we  know  that  there  exists  a 
constant  A/0  such  that  solutions  Ua.  Va,-  •  -.  Wh  are  bounded  by  A/„.  Define 

SJt)=  sup  A#0"-'max(||t/l«'>(-,s)  -  Uf'i-.  s)\\,*. 

0 

1 


(2.30) 


IIK,('”(*.s)  -  W-.*)lli-. 
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We  obtain  from  (2.29c) 

wr) 

(2  31)  =  "  mWm+1)  -  Wlm*U)  +  -  U^)V^wr  u 

+  mUil){Vaa)  -  Ki(1))W/jm  l) 

The  right-hand  side  of  (2.31),  denoted  by  h(x ,  t),  is  bounded  by 

PKOlL-  ^  w^m  +  i(0  +  2mA/0m St(/)  +  ^  4 m8m  +  1(r). 

By  the  assumption  of  Theorem  2.5,  lVjm>(x,0)  =  fV/jm\x,0).  Integration  of 

(2.31)  along  its  characteristic  line  will  give 

ll^r’KO  -  g4m/'Sm^(s)ds. 

Jo 

Moreover  one  can  show  that,  for  1  <;  p  g  m, 

-  W',(-'t)HL.£4’»fo'8m  +  i(s)ds. 

A  similar  result  applies  to  the  U  and  V  components.  Thus  we  conclude  that 

(2.32)  Sm(t)^4ntf'8m+l(s)ds. 

Jo 

Note  that  8m(t )  ^  A/0m.  By  induction  one  can  show  that 

(2.33)  6m(t)^4*^4l)T-HMo'”+‘  for  /c  =  0.1.--. 

For  T  <  1/4 M0,  the  right-hand  side  of  (2.33)  tends  to  zero  as  k  -*  oo.  Hence 
8m(t)  =  0.  This  proves  the  theorem  for  the  case  I. 

Case  II:  a  is  a  rational  number,  a  =  m/n.  We  assume  that  \m/n\  <  1.  Let  us 
choose  initial  data  of  type  a  to  be 

Ua{x,  ^,0)  =  sin(2?r(m  +  n)y)  +  1.0, 

Va(x,  y,0)  =  cos(2 ir(m  -  n)y)  +  1.0, 

Wa(x,  ^,0)  =  8.0(1  +  0.5  sin(4wn>’)). 
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and  choose  initial  data  of  type  b  to  be 

Ub(x,  y,(j)  =  sin(2w(m  +  n)y)  +  1.0, 

Vb(x ,  y,0)  =  cos(2w(m  -  n  )_y)  +  1 .0, 

Wb(x<  y,0)  =  8.0(1  +  0.5  cos(47rny )). 

Then  equalities  (2.28a-c)  are  satisfied  at  time  /  =  0.  However,  direct  computa¬ 
tion  shows  that 

f^Wa(x,  y,0)~J”ua(x,  y  +  (m/n  -  l)z,0) 

‘Va(x,  y  +  ( m/n  +  l)z.0)  d: j  dy  =  9.0. 

and 

f^  Wh(x,y.0)~fyh(x.y  +  (m/n  -  l)z.0) 

• Vh(x ,  y  +  (m/n  +  1)2.0)  dz |  dy  =  8.0. 

Thus  for  the  above  choices  of  initial  data,  we  get 

(2-34)  ^  J^(Wa(x,  y,t)2  -  Wh(x,  y\t)2)  dy  =  2.0  at  /  =  0. 

and  we  conclude  from  (2.34)  that 

f  (  Wa(x,  y,  t )2  -  Wh(x,  v,  t)2)  dy  >  0  for  t>  0  small. 

Jo 

Therefore  (2.28c)  cannot  be  true  for  all  t  >  0.  This  completes  the  proof  of 
Theorem  2.5. 

3.  Discrete  Boltzmann  Equations  with  Multiple  Scale  Initial  Data 

We  choose  the  simple  Carleman  model  to  illustrate  the  results.  Generalization 
to  more  complicated  models  follows  directly.  The  Carleman  equations  are  given 
as  follows  (see  [6]): 
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We  assume  that  the  initial  values  are  given  by 
(3.2)  u(*>0)  =  «0|x,  v(x,0)  = 


where  «0(x,  y,,  y2)  and  v0(x ,  y,,  y2)  are  1-periodic  functions  in  each  y,  variable. 

Case  1:  e2Ai  0  as  e2  -*  0.  The  homogenized  equations  for  (3.1 )— (3.2)  are 
given  by 


(3-3b)  TT  ~  Ti  +  yl  ~  OoU{x' y"  yi't)2  dy' dy 2  =  °' 

with  initial  data 

(3.4)  U(x ,  y,.  y2,0)  =  u0{x.  y,,  y2),  V{x.  y,.  y2.0)  =  v0(x.  y,.  y2). 

Case  2:  £i/e2  -*  a,  an  irrational  number,  and  |e,/c2  -  a|  g  c(Cj)r  with  r  >  1. 
Then  the  homogenized  equations  for  (3.1 )— (3.2)  are  given  by  (3.3)— (3.4).  the  same 
homogenized  equations  as  for  the  case  1. 


Case  3:  e,/e2  -*  m/n  *  0,  and  ^/e,  -  m/n\  s  c(E,)r  with  r  >  1.  Then  the 
homogenized  equations  for  (3.1 )— (3.2)  are  different  from  (3.3)— (3.4).  They  are 
given  by 

(3.5a)  +  +  V^x'  zA"i/n):,t)2  dz  =  0. 

(3.5b)  ~  "  It  +  vl-\iy(x<:Am/n)z,t)2dz  =  0. 

with  initial  data 

(3.6)  U(x ,  yx,  y2,0)  =  u0(x.  y,.  y2),  K(a.  y,,  y2,0)  =  d0(.v,  y,.  y2). 


Theorem  3.1.  Suppose  that  u0(.x,  v,.  y2),  i>0(x,  y,.  y2)  arc  bounded,  non¬ 
negative  and  continuously  differentiable.  Then  we  have 


u(x,t)  -  U\  x. 


0  as  Ej.  e2  -*  0, 


v(x,  t) 


X  +  t 
*1 


JC  +  /  \ 

— •') 


0  as  Ej ,  e2  — *  0. 


strongly  in  the  Lx-norm. 


488 


T.  Y.  HOU 


Remark  3.1.  The  above  homogenization  results  can  be  generalized  to  the 
case  where  the  initial  data  are  of  finitely  many  scales.  For  instances,  if  the  initial 
values  for  (3.1)  are  given  by 


«(x,0)  = 

o(x,0)  =  o0(x,~,^,^J, 


and  if  e1/e2  -*  m/n  *  0,  e3/e2  -»  0,  then  the  homogenized  equations  will  be 
fff  +  +  ~  n  f0  /0V(*-  Z1'(W/")21’  22'  t)2dz2dzl  =  0, 


dt 


JZ  +  yl  ~  i  /0  f^U(x,zu(m/n)zu  z2,  t)1  dz  2  dzx  =  0, 


with  initial  da?:. 


U{x,  y2,  y3, 0)  =  u0(x,  y \ ,  y2,  y3). 


>V  y2,  y^O)  =  v0(x,  >>„  y2,  y3). 


To  prove  Theorem  3.1,  we  need  a  few  technical  lemmas. 

Lemma  3.1.  Suppose  that  f{x,  y,  z)  is  continuously  differentiable  and  is  \-peri- 
odic  in  y  and  z.  Assume  that  t2/ei  0  as  Ej  -»  0;  then,  for  any  constants,  a,  b 
and  x0, 


(3.7) 


lim 

e\  “♦  0 


_j _ 

(b  -  a) 


~  j  -  fj(x'  y' z )  dydz] dx  =  °- 


Proof:  Arguing  as  in  Lemma  2.2,  we  need  only  to  prove  the  lemma  in  the 
case  when  /  is  independent  of  the  first  variable.  We  assume  for  simplicity  that  / 
has  been  normalized  so  that 


(3.8) 


ff/iy '  z)dydz  =  0. 
Jo  Jo 


HOMOGENIZATION  FOR  SEMIL1NEAR  HYPERBOLIC  SYSTEMS 


489 


For  simplicity,  we  let  x0  =  0.  By  change  of  variable,  the  left-hand  side  of  (3.7) 
becomes 


Since  f(y,z)  is  assumed  to  be  1-periodic  in  y  and  z,  the  first  term  on  the 
right-hand  side  of  (3.9)  is  bounded  by  ce,  in  light  of  (3.8).  Define 

g(y>  2)  sf(y>  z)  -  pf(y,z)dz. 

Jo 

Then  g(y, z)  is  1-periodic  in  z  and  satisfies  /0’g(>\  z)  dz  =  0.  The  last  term  in 
(3.9)  is  then  bounded  by 


(3.10) 


—  f  /e'g\  X,—)dX 

~  °Ja A,  '  £2  ) 


g  sup 


/ 


n+  1 


g\ 


dx 


+  0(e. 


Applying  Lemma  2.1  to  J„n+1g(x,  x/le^/e^)  dx,  we  conclude  that  the  left-hand 
side  of  (3.9)  converges  to  zero  as  -»  0.  This  completes  the  proof  of  Lemma  3.1. 


Lemma  3.2.  Suppose  that  f(x.  y,  z)  is  continuously  differentiable  and  is  l -peri¬ 
odic  in  y  and  z.  Assume  that  et/e2  =  a  4-  0((e,)r)  with  a  *  0  and  r  >  1;  then, 
for  any  constants  a,  b  and  x0, 


dx 


(3.11) 


T^fffahl(f{x'z'{m/n)zS>dzdx 


if  a  is  irrational, 
if  a  =  m/n. 


Proof:  Arguing  as  in  Lemma  2.2,  we  need  only  to  prove  the  lemma  in  the 
case  when  /  is  independent  of  the  first  variable.  For  simplicity,  we  assume 
x0  =  0.  By  change  of  variable,  we  obtain 
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The  second  term  on  the  right-hand  side  of  (3.12)  is  bounded  by 


(3.13) 


fl>  Ai 

*1 

/  \y\ 

- - a 

p 

■'a/e, 

fc2 

dy  ^  Ce[ 


r-  1 


For  the  first  term  on  the  right-hand  side  of  (3.12),  if  a  is  irrational.  Lemma  2.2 
implies 


(3.14)  lim  ——  (h/t'f{y,ay)  dy  =  f  ff(y,  z)  dydz. 

*,-O0  QJa/e,  Jo  Jo 


If  a  =  m/n,  then  f(y,(m/n)y)  is  an  n-periodic  function  in  y.  Thus  we  have 
<315) 


Lemma  3.2  then  follows  from  (3.10)— (3.15). 


Proof  of  Theorem  3.1:  The  solutions  of  the  Carleman  equations  are  known 
to  be  bounded  for  all  time  for  bounded  non-negative  initial  data  (see  [10],  [13]). 
Thus  Theorem  3.1  can  be  proved  by  using  the  similar  techniques  we  use  in  the 
proof  of  Theorem  2.1  and  the  two  lemmas  above.  We  omit  the  proof. 

Remark  3.1.  If  el/e2  -*  a  *  0  and  the  rate  of  convergence  is  of  order 
0((ei)r)  with  r  ^  1,  then  Lemma  3.2  may  not  hold  in  general. 

For  examples,  suppose  that  e,/e2  =  1  +  Ej  (r  =  1).  We  take  /( y,  z)  = 
cos(27ry)sin(2wz),  x0  =  0,  a  =  0,  b  =  Then  a  direct  calculation  shows  that 

2 /1/2/f  7-,  f )  dx  ~  dz  =  1/”'  +  °^Ei)  *  °- 

J0  \£i  e2l  Jo 


4.  General  Model  for  a  Gas  with  Discrete  Velocity  Distribution 

Consider  a  gas  composed  of  identical  particles  of  mass  m.  The  velocities  of 
these  particles  are  restricted  to  a  given  finite  set  of  p  vectors:  Uj,  ,  up. 
Nt  =  TV  (x,  t)  denotes  the  number  density  of  particles  with  velocity  u,  at  the  point 
x  and  at  the  time  t. 

We  consider  the  binary  collisions  only.  Denote  by  u,  and  u,  the  velocities  of 
two  molecules  before  an  encounter,  after  the  encounter  these  molecules  have 
velocities  and  u,.  They  must  satisfy  the  following  two  relations  expressing  the 
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conservation  of  momentum  and  the  conservation  of  energy: 

U,  +  u,  =  u*  + 

M,|2  +  |u,|2  =  |U*j2  +  |U,|2. 


A  transition  probability  A*J  is  associated  with  the  collision  “u,,  u,  to  u*.  u/’. 
A^NNj  is  the  number  of  collisions  u1(  to  u*.  u,  per  unit  time  and  unit 
volume,  satisfying  the  particle  indistinguishability: 


tki 


Then  the  Boltzmann  equation  is  replaced  by  a  system  of  p  nonlinear  partial 
differential  equations  (see  [10]): 

(4.1)  ^+u,.VjV(  =  \  £  Ay).  ;  =  \,2,  -  ■  ■ ,  p. 

J.k.l 

We  are  especially  interested  in  the  case  when  the  initial  values  are  of  the 
form: 

(4.2)  Ni  0(x)  =  J),  i- 1, 2. •••,/>. 

where  <p,(x,  y)  are  1-periodic  functions  in  each  component  of  y. 

Assumption.  Let  { y,,  1  s/sn)  be  a  sequence  of  non-zero  numbers. 
Suppose  that  { y x , -  -  - ,  yA }  is  the  largest  linearly  independent  set  of  {YiC  •  Y«} 

among  the  integers  in  the  following  sense: 

We  say  that  { Y;, ■  -  • ,  y^ }  is  linearly  independent  among  integers  provided 
m  *  (Yi>-  •  •»  Y*)  =  0  f°r  some  integer- valued  vector  m  implies  that  m  =  0. 

Since  { YiC  •  •,  Y* }  is  the  largest  linearly  independent  set  of  { Yi,  ■  ■  • ,  y„ },  there 
exist  integers  /  and  m  such  that 

i  k 

(4.3)  Y,  =  —  L  t,,jYj  for  i  =  k  +  1. •  ••,«. 

'  j- 1 

We  denote  by  Nk  the  smallest  common  integral  multiple  of  mk  +  l,-  ■  •,  m„. 

Lemma  4.1.  Suppose  f(xv  •  •  •,  xn)  is  a  continuous  function ,  and  is  l-periodic 
in  each  component.  Let  { y,,  1  S  i  £  «}  be  defined  as  in  the  above  assumption. 
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Define  function  g(z2,-  ■  zk)  by 


g{z2r  ■  •-  zk)  =  fM/\  x2  +  Y22i  +  z2 * '  ■••**  +  Y*’i  +  rA. 

‘'n  \ 


*  / 

E*  *  +  1  ,  ffl 

— - z- 


2 


m 


*+1 


■■•.*»  +  Y«*l  +  L  JZrZm)dz  I' 

m=2  n  / 

Then  g(z2,-  ■  ■,  zk)  is  a  continuous  ^-periodic  function  in  each  zm  for  2  g  w  <;  A:, 
and 


[T/MJ-l 

//=y  I  «(7>2A/.”-.yY*^)  +  0(l/r). 
y-o 

Suppose  that  there  exists  an  integer-valued  vector  m  e  Z*  such  that 
^2(^2^)  +  •  ••  +"1*:(y*A/)  =  -m,ez. 

Recall  that  M  =  A^/y,.  So  we  have 

m,Yi  +  m2Nk  y2  +  ■■■  +  mkNkyk  =  0. 

The  assumption  on  {  y,;  1  g  i  ^  k }  implies  that  m,  =  0  for  1  g  /  g  k.  Thus 
we  conclude  that  {y,M\  2  ^  ^  k)  are  linearly  independent  among  integers  in 
the  sense  of  (1).  Applying  the  theorem  of  “Ergodic  Translation  of  Tori”  in  [1]  to 
g(z2,-  •  •,  zk)  with  u  =  (y2M%-  ■  -,ykM),  we  obtain 

(4.8)  lim  //=  — y-+—  fNk  -..  fNkg(z2,---,zk)  dz2  •••  dzk. 

T—  oo  (Nk)  J0  J0 

Thus  we  prove  (4.6)  by  expressing  the  right-hand  side  of  (4.8)  in  terms  of  /.  This 
completes  the  proof  of  the  lemma. 

Now  we  can  define  the  homogenized  equations  for  the  general  discrete 
Boltzmann  equations. 

Suppose  is  a  continuous  1-periodic  function  in  To  1  S  «  £  «■ 

Let  Y  =  (Yi>‘  ’  *>  Y„)- 

Definition.  Define  functional  F[f;  y]  as  the  limiting  function  on  the  right- 
hand  side  of  (4.4),  (4.5)  and  (4.6)  corresponding  to  the  cases  k  =  1,  k  =  n  or 
1  <  k  <  n,  respectively. 
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For  x  and  t  fixed,  we  regard  A/y(x,  y,  /)  as  a  function  of  y  alone.  Define  the 
functional  G,(My,  Mk)  as  follows: 


y.  /),  A/t(x,y,  /))  =  { 


{MjMk 

M/lM.A u,-u*)] 
M,F[My;(u,-uy)] 

(u,  -  uy.u,  -  uj] 


if  J  =  *  =  ». 
if  /  =  i,  k  *  /, 
if  j  *  i,  k  =  i. 
if  j  ¥=  i,  k  *  i . 


where  F[MjMk\  (u,  -  uy.  u,  -  u^)]  corresponds  to  the  limiting  function  in 
Lemma  4.1  with  /  given  by 

/( y  +  +  pr)  =  Mj(’.y  +  (u,  -  u ,)/.  ‘)Mk(  *.z  +  (u,  -  uA )/,  •). 

Then  the  homogenized  equations  are  given  by 
S  M  1 

(4.9)  -jf  +  u,- VXM,  =  ^  M#)  -  M,)). 

J.  k.l 

with  initial  values 

(4.10)  M,(x.y.O)  =  </>,(x.y),  '  =  1 . •  •  ■ .  ^ ■ 


Theorem  4.1.  Let  M,  be  the  solutions  of  the  homogenized  equations 
(4.9)-(4.10).  Then  we  have 

E  k(x,()-W,  x, - - — -»  0  as  e  *  0, 

1-1  II  '  Mlz«(R!;[0.r|> 

provided  that  equations  (4.1  )-(4.2)  and  (4.9)— (4.10)  have  bounded  solutions  for 
OS^T. 

Proof  of  Theorem  4.1:  Using  Lemma  4.1,  we  can  prove  Theorem  4.1  in  a 
similar  way  to  Theorem  2.1.  We  omit  the  proof. 

Remark  4.1.  The  global  existence  theory  of  equations  (4.1)— (4.2)  has  been 
investigated  by  several  authors  (see  e.g.  [12],  [10),  [5],  [13]  and  [2]).  The  general 
results  in  [10]  indicate  that  the  bounded  global  solutions  exist  for  long  time  if  the 
initial  data  are  small  in  some  sense.  The  local  existence  of  the  homogenized 
equations  (4.9)  can  be  obtained  by  classical  analysis.  By  combining  the  known 
global  existence  results  for  (4.1)  with  Theorem  4.1,  we  can  show  that  bounded 
global  solutions  of  (4.9)  exist  as  long  as  the  global  solutions  of  (4.1)  exist. 

Remark  4.2.  The  techniques  we  use  here  in  deriving  the  homogenized 
equations  can  be  used  for  more  complicated  models.  For  examples,  our  tech- 
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niques  will  apply  to  the  semi-linear  hyperbolic  systems  of  the  form 
dN 

-jf  +uI-V*I«/I(Af„- •■,**),  /  = 

with  oscillatory  initial  data 

^,(x.O)  =  Af,.0(x,  j),  /  =  1.2.  -.k. 

if  the  u,  are  constant  vectors  and  the  /,  are  smooth  and  separable  functions  in 
Nrj  =  h---k. 
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